manchester 
particle physics 



MC-TH-97/04-rev 
MAN/HEP/97/3-rev 
June 11, 1997 



Estimating the effects of Bose-Einstein correlations 
on the W mass measurement at LEP2 

V. KartvelishviliQ 

Department of Physics and Astronomy, University of Manchester, UK, 

R. Kvatadze 1 
Institut de Physique Nucleaire de Lyon, France 

and 
R. M0ller 

Niels Bohr Institute, University of Copenhagen, Denmark. 



Abstract 

The influence of Bose-Einstein correlations on the determination of the mass of 
the W boson in e + e~ — ► W + W~ — ► 4 jet events at LEP2 energies is studied, using 
a global event weighting method. We find that it is possible to keep the systematic 
error on the W mass from this source below 20 MeV, if suitable precautions are 
taken in the experimental analysis. 
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1 Introduction 



The accurate determination of the mass of the W boson is expected to be one of the most 
important standard physics results to be obtained from LEP2 measurements. Since M^y 
is one of the key parameters of the electroweak theory, precise knowledge of the W mass 
allows one to test the Standard Model (SM), and together with the top quark mass can 
be used to constrain the allowed range of the Higgs boson mass in the SM, or restrict the 
parameter space of other "new physics" . The present value of Mm is based on the direct 
measurements in pp interactions at CERN |l] and at the Tevatron [J. The combined 
result from these studies is Mw = 80.33 ± 0.15 GeV @]. Data from LEP2 recorded in 
1996 at 161 and 172 GeV are expected to lead already to a comparable error when they 
are fully analysed and the results of the four LEP experiments are combined. When all 
existing data from the CDF and DO experiments (more than 100 pb _1 per experiment) are 
analyzed, it is envisaged that the error on the W mass will be reduced to about 70 MeV. 
On the other hand, an indirect determination of Mm from a SM fit using LEP and SLC 
measurements give an error in the W mass of approximately 40 MeV (for a Higgs mass 
of 300 GeV) |J, thus setting the scale for a significant new test of the model. 

Three different methods have been proposed for the determination of the W mass in 
e+e" annihilation at LEP2 0, ft g: 

- The threshold cross-section measurement of the process e + e~— > W + W _ . It has been 
shown that maximum sensitivity to M^y is achieved for the energy ^/s ~ 161 GeV. 
Now these measurements are complete and the error in the W mass is about 450 MeV 
for a single experiment ||, with a combined result 80.4 ± 0.22 GeV |l( 



- The measurement of the charged lepton end-point energy. This gives an estimated 



error exceeding 300 MeV [pi 11] for the total integrated luminosity expected at LEP2 



and is therefore not competitive with the other two methods. 

- The direct reconstruction of the Myy from the final state particles. This is consid- 
ered to be the most promising method. It relies on the use of the four constraints 
from the energy- momentum conservation. An additional constraint can be added 
by the assumption that the two W bosons have equal masses. Two decay channels 
- qqqq (four jets) and qqlv (two jets plus leptons) — can be used in this analysis. 
The estimated total error for the two channels combined is 34 MeV, assuming four 
experiments each collecting 500 pb _1 data at ^fs = 175 GeV ||. 

For the process e + e~— > W + W~ at LEP2 energies, the typical separation of the two 
decay vertices of Ws in space and time is of order of 0.1 fm, much smaller than the 
hadronization scale (~ 0.5 fm). Thus, when both Ws decay hadronically, the hadroniza- 
tion regions of the W + and W~ overlap, and colour interconnection |i~2|, |I"3f during the 



hadronization as well as Bose-Einstein (BE) correlations [14] between identical bosons 
among the decay products of the two Ws can couple the two systems and thereby affect 
the W mass measurement. 

At present it is not excluded that these effects could each contribute about 100 MeV to 
the theoretical uncertainty, which would make the four jet channel essentially useless for 
the W mass measurement. Hence, apart from the intrinsic interest in these phenomena 
the studies of colour interconnection and Bose-Einstein correlations are very important, 
and reliable estimates of the size of their effect on the reconstructed W mass are highly 
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desirable in order to evaluate their contributions to the systematic error in this measure- 
ment. 

The first attempt to investigate the influence of Bose-Einstein correlations in the de- 



termination of W mass was made in [14]], where the LUBOEI algorithm, implemented 



in the JETSET Monte-Carlo program ||15|| , was used. The basic assumptions of this ap- 
proach are that BE effects are local in phase space and do not alter such characteristics 
as the event multiplicity or the cross section. The momenta of the bosons produced at 
the hadronization stage are shifted by amounts calculated to reproduce the two particle 
correlation function expected for a source with a gaussian space-time distribution. This 
procedure, however, does not preserve momentum and energy simultaneously. The vi- 
olation is not large, however, and energy conservation is restored in an ad hoc way by 
rescaling all particle momenta with a common factor. The latter procedure introduces 
an artificial shift in the W mass, even if there are no Bose-Einstein correlations between 
the particles coming from different Ws. The procedure also tends to make the jets more 
narrow, thereby decreasing the individual jet masses. Assuming a source radius of 0.5 
fm, the corrected mass shift was found to be 95 MeV at ^/s = 170 GeV [13]. The shift 
increases with increasing energy and decreasing source radius. 

In the following analysis we will address the problem of Bose-Einstein correlations, 
but using another approach, which is based on assigning weights to the simulated events 
according to the momentum distributions of final state bosons. In the global event weight 
scheme a shift in the W mass can arise, if the event weight depends on Mm. The 
use of such global event weights is not straightforward and was discarded in the above 
study on fairly general grounds. Our aim is to reinvestigate the suitability of the method 
and to estimate approximately the systematic error on the W mass from Bose-Einstein 
correlations. The use of global event weights is complementary to the local reweighting 
scheme of ref. jL4j in the sense that here the kinematical properties of the events are 
preserved, while all probabilities and multiplicities may change. We have not found a 
unique solution to the problem of assigning the weights, however, and have used a number 
of different weighting schemes to assess the range of effects that can be expected from 
Bose-Einstein correlations. 

We have used the PYTHIA event generator |15] to generate Monte Carlo samples 
of e + e~ — > W + W~ — > 4 jet events. A basic assumption here is that hadronic W and Z° 
boson decays are sufficiently similar, so that by using the tuning of the Monte-Carlo model 
parameters that reproduces the experimental data from Z° decays at LEP, Bose-Einstein 
effects in single W decays are already effectively taken into account in properties such as 
multiplicities and single particle momentum spectra. Hence, only the correlation between 
identical particles from different Ws have been included in the weight calculation. In 
order to check the self-consistency and inherent systematic errors of the method we have 
applied the same weighting method to the well studied process of Z° hadronic decays. 

Since we want to estimate the maximal systematic error on the measured W mass due 
to the mass shift arising from Bose-Einstein correlations, we have not taken into account 
smearing due to experimental resolution and acceptance, reconstruction method etc. 

Various possibilities of constructing event weights are discussed in the following sec- 
tion. Some consequences of these weighting schemes on certain measurable quantities 
at the Z° peak are considered in Section 3, together with possible constraints on BE ef- 
fects from Z° physics. In Section 4 the influence of BE effects on the process of W pair 
production at LEP2 is analyzed. 
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2 Event weighting schemes for Bose-Einstein effects 



The Bose-Einstein effect corresponds to an enhancement in the production probability of 
identical bosons to be emitted with small relative momenta, as compared to non-identical 
particles under otherwise similar conditions. Experimental data are usually analysed in 
terms of the correlation function defined as the ratio of the two-particle probability density 
to the product of the corresponding single particle quantities: 

where pi is the four- momentum of particle i. Assuming a spherical space-time distribution 
of the particle source, the correlation function takes the form: 

C{Q) = 1 + Xp(Q) (2) 

where Q is the four-momentum difference, Q 2 = — (pi — P2) 2 , and p is equal to the 
absolute square of the Fourier transform of the particle emitting source density, with the 
normalization condition p(0) = 1. The parameter A, known as the incoherence parameter, 
takes into account the fact that for various reasons, the strength of the correlations can 
be reduced. 

Often a gaussian model is assumed for the source density, which leads to 

p(Q) = exp(- J R 2 g 2 ) (3) 

where R is the radius parameter. Experimentally, this gives a reasonable description of 
Bose-Einstein correlations in many types of collisions. In e + e~ collisions typical values of 
R and A are respectively 0.5 fm and 0.3 (1.0 for directly produced pions) (T^, fH| . 

In order to simulate the effects of Bose-Einstein correlations, we have in this study 
chosen to use the event weighting method. The method arises very naturally in a quantum 
mechanical approach, where the weight can be constructed as the ratio of the square of 
the symmetrized multiparticle amplitude to the square of the non-symmetrized amplitude 
corresponding to the emission of distinguishable particles. The use of global event weights 
leads to a number of conceptual and computational difficulties, however, which must be 
overcome for any quantitative conclusions to be drawn and which we will attempt to 
address below. 

There are several possibilities to construct event weights. One way of forming the 
weight is to take a product of enhancements C(Q) for all pairs of identical bosons in the 



event WE\: 



Vi=U C(Q hi2 ) (4) 

For high multiplicity events this weight can become extremely large, so that a few 
such events dominate the weighted distributions and lead to non-realistic distributions. 
The event weights therefore have to be regularized in some way. In order to keep the 
statistical error at a reasonably low level, we have chosen to discard events with very high 
weights (higher than some V max ). This, however, has the unpleasant effect of making the 
results Vm^-dependent. We have tried to circumvent this difficulty by analyzing the V max 
dependence and extrapolating the results to V max — > 00. 
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Figure 1: Distributions of the weights V\ (dashed line), V 3 (full line) and V 5 (dotted line) 
for WW events at 175 GeV. 

One can also rescale the weight of the event using a single constant Wq: 

V 2 = V x /w% (5) 

where wq is a constant slightly larger than 1, and n is the number of pairs in the event (i.e. 
the number of terms in the product in ([|)). The value of wq is chosen to keep the average 
multiplicity reasonably close to its value before event weighting P0[| . For a constant V max , 
we have found, however, that the method gives rise to numerical difficulties, stemming 
from the fact that increasing u> brings in more events from the high weight tail of V\, 
which leads to large fluctuations in the multiplicity and the average event weight. Our 
results for the shifts in multiplicity and M^y using V 2 are roughly consistent with those 
found using Vj., and V 2 will not be discussed further. 

A problem with these weighting methods is that, like the local reweighting scheme 
of ref ||14|| , they emphasize narrow jets and thereby introduce artificial correlations also 
between non-identical particles. In order to counteract this, one can use the weight 
calculated with non-identical pairs to rescale (01): 

V 3 = V x fVZ /m (6) 

where Vq is the weight calculated according to but for non-identical bosons in the same 
event, while n and m are the numbers of identical and non-identical pairs, respectively. 
This also leads to a better numerical behaviour, as illustrated in fig. 1, which shows the 
distributions of V\ and V3 for simulated WW events at \fs = 175 GeV. The high weight 
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tail is much less pronounced for V3 than for V\. Both fall off to a good approximation 
as inverse powers, with exponents —2.6 and —1.4, respectively, which makes it plausible 
that the sum of all weights converges. 

A different method of constructing the event weight, which is closer to a full quantum 
mechanical treatment, starts from the introduction of a symmetric amplitude, which has 
n! terms [pOfl . This leads to a weight: 



v 4 = A fc /2 p (Q ln )p(g 2l2 )---p(g m j, (7) 

permutations 

where k is the number of times when the first and second indices differ. For Q — 0, 
A = 1 and n identical particles, equation (B) gives a weig ht of 2"( n ~ 1 )/ 2 , while eq. © 
results in the correct value n\. However, for typical hadronic configurations this difference 
is much smaller, and (|7|) is rarely used because of computational difficulties. We had to 
restrict ourselves to events containing no more than 8 identical particles of each kind. 
This limitation is too restrictive already at Z° energies, where the number of identical 
boson combinations is lower than in W + W~ events, and rejected about 50 % of events in 
W pair production, making it effectively useless for the latter case. We used it only to 
check that it gave essentially consistent results with our other methods for events with 
low multiplicity. 

The event weights defined above were based on the gaussian parametrization of the 
particle emitting source. This implies that V\ > 1 for all values of Q. 

This is not true in all models, however. In addition to the above weights, we have 



therefore also studied a different pair weight, inspired by the weight used in where p 
in (g) is not required to be always positive: 

cos(£Qi?) 

piQ) = ^ (8) 

For £ close to 1, this is very close to (|3|) apart from becoming slightly negative at large 
Q. The corresponding weight V5 is built in analogy with but with the gaussian @ 
replaced by (|8|) (the dotted line in Figure 1). We find that £ = 1.15 leads to a good overall 
description. Due to the better numerical behaviour of this weight function (the exponent 
in the power fit is —2.4), we were able to apply V5 without further rescaling. 



3 Influence of event weighting on Z properties 

Various measurable properties of the Z° will be affected to different extent, if one intro- 
duces event weights into the simulation of its hadronic decays. Since the partonic states 
before hadronization are known to be well described by perturbative calculations, which 
do not take into account Bose-Einstein correlations, uncritical application of event weights 
may lead to large inconsistencies with e.g. measured branching ratios and relative frequen- 
cies of jet multiplicities etc. In order to see how serious these effects are and to judge 
what consequences this has for the analysis of the WW events, the precise experimental 
data from Z° decays can be used to check the event weighting schemes of Bose-Einstein 
correlations for W pair production at LEP2. We have simulated 3 x 100000 hadronic 
events at y'i = M% and M% ± 2 GeV. Table 1 presents the differences for charged parti- 
cle multiplicity, shift of Z° peak position in hadronic vs leptonic decay modes, branching 
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v 1 


v 3 


\ 




A(n ch ) 


3.7 ±0.5 


1.3 ±0.2 


1.8 ±0 


2 


AM Zo , MeV 


8±3 


0±3 


1 ± 


4 


AR C , % 


-3 ±2 


-2 ±2 


± 


2 


AR b , % 


-26 ±3 


-11 ±2 


-5± 


2 


A 3jet/2jet, % 


80 ±20 


20 ±5 


20 ± 


5 



Table 1: Differences in charged multiplicity, peak mass of Z°, branching fractions 
and three-to-two jet event ratio, between weighted and non-weighted events, for various 
weighting systems described in the text. 

fractions for charm and beauty decays (R c and Rb) and the ratio of three to two jet events 
with and without event weighting. 

This analysis resulted in the following: 

- The average charged multiplicity has changed. The weight V%, which was not 
rescaled, leads to the largest increase when compared to the unweighted results, 
while both V3 and V5 give a smaller increase around 1.5. In all these cases, the 
change can be accommodated by retuning of the parameters in the simulating pro- 
gram. 

- In principle, event weighting can result in a shift of the Z mass peak. However, only 
Vi yielded a shift of a few MeV, while for V3 and V5 the shift is essentially zero. We 
have not found any significant change of the Z° width. 

- The pattern of heavy and light quark fragmentation is rather different. Heavy quarks 
produce significantly less pairs with small Q, and all BE effects in this approach are 
less pronounced for heavy quarks. Heavy quark events thus obtain smaller average 
weights, which result in changes shown in Table 1. Note that the effect for c-quarks 
is diluted because fe-quark events reduce the overall average weight. Further study of 
this effect lies beyond the scope of this work. In order to exclude this artificial flavour 
dependence in W decays, the weighting and rescaling was performed separately for 
the different decay modes of the Ws. 

- The weighting resulted in a substantial increase of jet activity, as measured by the 
three to two jet event ratio. This is however difficult to quantify because of its 
dependence upon the jet finding algorithm and its parameters. The numbers shown 
in Table 1 were obtained using LUCLUS with default parameters, corresponding to 
fairly narrow jets. The effect decreases for broader jets and in any case is much less 
pronounced in WW production, so we did not attempt to correct for it. 

We have considered maximum BE correlations, A = 1, but only for pions and kaons 
originating from the sources with decay lengths cr < 10 fm, and A = otherwise. The 
source radius R was taken to be equal to 0.5 fm everywhere. In Z° decays at LEP1 one 
observes A ~ 0.3 if all particles are considered, 0.4 if only pions are taken into account 



and 1.0 for directly produced pions [IS, IS], and R « 0.5 fm 
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Figure 2: Reproduced correlation functions in Z° events using the weights V\, V3 and V5. 
The dashed lines show the result of a fit to real data from hadronic Z° decays JL9 . 



The reproduced correlation functions for the three weighting schemes, V\, V3 and V5 
are shown in figure 2. We have divided the ratio of the weighted same-sign particle 
distribution to the unweighted one by the similar ratio for opposite-sign particles. This 
is equivalent to using the Q-distribution of opposite-sign pairs as a reference sample and 
correcting for effects of particle selection and resonances by dividing by the same ratio 
in simulated events without BE correlations — a common procedure in experimental 
analyses. Also shown are fits to the form 



N(l + PQ)(1 + Xexp(-Q 2 R 2 )) 



(9) 



which is often used to parametrize the experimentally observed correlation function in Z° 



decays [18, 19]. The resulting values of the parameters for a fit range of 0-2 GeV in Q 
are shown in Table 2. The dashed line in the figure represents the result of a fit to the 
correlation function of all particles observed in real data from hadronic Z° decays [|nj . 

V\ does not reproduce the correlation function well, as it gives too small values for 
both the incoherence parameter A and the radius R. The other weight schemes give 
very reasonable descriptions resulting in input and output parameter values which are 
reasonably close to each other. We take the spread of the results using the different 
weighting schemes to be indicative of the systematic errors inherent in the method. 

Hence we conclude that, provided that the different quark final states (and possibly 
the final states with different number of jets) are treated separately, application of the 
global event weighting technique with rescaling of the weight (V3) or using the form (§) is 
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Weight 


Vi v 3 v 5 


A 


0.072 ±.004 0.292 ±.005 0.325 ± .005 


R (fm) 


0.387 ±.018 0.465 ±.005 0.455 ± .005 


(3 (GeV" 1 ) 


-0.029 ±.003 0.047 ±.003 0.080 ± .003 


N 


1.010 ±.004 0.925 ±.004 0.899 ± .004 



Table 2: The fitted values of the correlation function parameters A, R, (3 and N at the 
Z° peak. 



Weight 


Vy V 3 V 5 


A 


0.032 ±.004 0.101 ±.002 0.146 ±.003 


R (fm) 


0.322 ±.024 0.459 ±.009 0.457 ± .006 


(3 (GeV" 1 ) 


0.002 ±.003 0.013 ±.001 0.036 ± .001 


N 


0.997 ±.003 0.977 ±.003 0.949 ± .003 



Table 3: The fitted values of the correlation function parameters A, R, (3 and N for W 
pair production at <Js = 175 GeV. 

not inconsistent with LEP1 data at the Z°, whereas the direct application of the product 
of pair weights (Vi) should be treated with more care. 



4 W pair production 

PYTHIA 5.7 was used to simulate the process e + e _ — > W + W~ — > 4 jets, and the weighting 
schemes described above were applied to generate BE effects. 

Only correlations between identical bosons originating from different Ws were in- 
cluded, since BE correlations within a single W cannot lead to any change in M^y com- 
pared to the semileptonic decays e + e _ — > W + W~ — > 2 jets+lu. 

The first thing to study is the correlation function and the compatibility of input and 
output values indicating the self-consistence of the method. Table 3 refers to 175 GeV, 
and contains the results of the fits to the correlation functions for each weight used (as 
above for Z°) using the cut off V max = 80. Numbers for 192 GeV are very similar. The 
errors given in the table are statistical only. From the variation of the fit results with 
Vmax, we estimate that the systematic error is 0.05 on A and 0.05 fm on R. 

One sees that for V\, R is somewhat lower that the input value, while V 3 and V 5 give 
values quite close to 0.5. It is worth noting that the values of A are significantly smaller 
than at the Z peak, essentially because in the WW case we have included only correlations 
between pairs from different Ws. 

Next, the information from the Monte-Carlo was used to assign each final particle to 



the W + or the W~, as in [Q. Ws with the mass values in the interval 70 GeV < < 
90 GeV were studied at 175 and 192 GeV to assess the energy dependence. At each energy, 
10 5 events were generated, which is about an order of magnitude higher than the expected 
statistics of all four LEP experiments combined at 500 pb _1 integrated luminosity per 
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Vi 


v 3 


v 5 










175 GeV 


3.8 ± 0.5 


1.8 ± 0.2 


1.0 ± 0.2 


192 GeV 


3.7 ±0.5 


1.7±0.2 


0.6 ±0.2 


AM^ ean (MeV) 








175 GeV 


75 ± 15 


22 ± 11 


20 ± 14 


192 GeV 


92 ± 16 


34 ± 11 


38 ± 14 


AM^ (MeV) 








175 GeV 


12 ±9 


11 ±7 


4± 12 


192 GeV 


15 ±8 


13 ±7 


6±9 



Table 4: Values of differences in multiplicity and mass of the W boson for events with 
and without interconnecting Bose-Einstein correlations between the two Ws. 

experiment. In general, one expects that BE-induced effects in WW production should 
die out at high energies, as the overlap between the two W decay volumes decreases. This 
requires much higher energies than will become available at LEP2, however, and it is 
likely that the effect will increase with energy in the LEP2 range ||14|| . 

The mass distribution of W bosons was built with and without event weighting for each 
of the weights used, and the differences were calculated in the average charged multiplicity 
n c h, the mean W mass, M™ an , averaged over the whole interval 70 GeV < Mw < 90 GeV, 
and a fitted M$. The fit was performed using a relativistic Breit-Wigner shape with an 
s-dependent width, in the interval 80.25 ± 5 GeV, with 5 = 2 GeV. The results are 
presented in Table 4. We also performed the fit in wider intervals corresponding to 5 = 
4,6,8 and 10 GeV and observed a small additional increase in the mass shift of a few MeV, 
which saturated at the larger interval sizes, consistent with the fact that the tails of the 
Breit-Wigner distribution contain very little information about the value of the W mass. 

As mentioned above, for computational reasons, we were forced to cut away a tail of 
events with very large weights. We have tried to eliminate the dependence on the cutoff 
value, V max , by calculating the multiplicity and mass shifts for three values of V max (20, 40 
and 80) and then extrapolating to infinite cutoff. This method seems to be more reliable 
and less vulnerable to fluctuations than direct calculation with very high V max . Figure 3 
shows the values of the mean W mass, M™ an , and the fitted as functions of 1/V max - 
There is no indication in our investigations that the inclusion of the events with very large 
weights would change the estimated mass shifts by any significant amounts. For Vi, the 
extrapolated value depends on the specific way the extrapolation is performed. We have 
included this ambiguity into the error shown in Table 4. 

From these numbers we draw the following conclusions: 

- There is a clear correlation between the BE-induced shifts in the W mass and in 
the charged particle multiplicity in WW production: the larger is the increase in 
charged multiplicity, the larger mass shifts are expected. 

- Both V3 and V5 weights result in fairly small shifts, while still maintaining a good 
reproduction of the correlation function. They are well-behaved numerically and in 
our opinion give quite reliable estimates of the effect. However, we conservatively 
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Figure 3: Shifts in the mean and fitted W mass as functions of the inverse weight cut 
off, l/V max , for weight schemes V\, V 3 and V 5 , and y/s =175 GeV. 



take the spread of values using all three weighting schemes as indicative of the 
systematic errors inherent in our approach. 

The fitted value for the W mass is less sensitive to BE effects than the mean over the 
full distribution, which has been used to estimate the effect in previous investigations 
|T4]j . Our estimated values for the shift in the fitted mass are less than 20 MeV, 



implying that BE correlations are not too dangerous for the W mass measurements 
at the expected level of accuracy at LEP2. For the shifts in the mean W mass, we 
confirm the conclusions of fl4| and find values of the same general magnitude of a 



few tens of MeV, under similar conditions as investigated there. In all cases the 



shift is towards larger masses, as expected on general grounds |BL 14 



- For all weighting schemes, the shift in Mm increases with energy in the energy 
range considered, but the increase is fairly small. 

It is interesting to compare our results to the predictions based on the implementation 
of Bose-Einstein effects by shifting the momenta of final state particles [pi}] . There are 
several differences in the predictions of these two schemes. 

The most important difference is in the particle multiplicity: our approach naturally 
leads to an increase of the average number of particles due to Bose-Einstein correlations, 
while the momentum- shifting method assumes that the multiplicity is unchanged. Ex- 
perimentally, it is not known yet to what extent Bose-Einstein correlations might modify 
the particle multiplicities at high energies. 
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The energy dependence of W mass shift is different. The strong energy dependence in 
momentum-shifting scheme is a combination of two effects: the increase of the systematic 
shift for low momentum particles in the direction of smaller W momenta, and the differ- 



ences in momentum spectra of W decay products for various energies, as stressed in |14|. 
This seems to be less pronounced in the present approach. 

The present study confirms that the systematic effect of BE correlations on the W 



mass determination can potentially be quite large, as found in [13], although the actual 
values of the mass shift found here are somewhat smaller. The size of the shift is however 
quite sensitive to the procedure used to extract the value of the W mass. In particular 
we observe that a fit to the lineshape of the W mass distribution has a much smaller 
systematic error from Bose-Einstein correlations than the average mass, due to the fact 
that the main effect on M^y in our scheme arises from the tails of the mass distribution, 
which contain very little information about the peak position. Hence it does seem possible 
to keep the systematic error from this source below about 20 MeV. Careful work linked to 
the actual fitting procedures used by the LEP experiments is obviously needed in order 
to assess this in the individual cases and to optimize the analysis procedures. Since the 
value of the mass shift is always positive (as also expected on general grounds), a further 
reduction of the systematic error by a factor two is in principle possible by assigning the 
expected shift correction to Mw. 

The estimates in this paper were all made assuming full Bose-Einstein correlation 
strength (A = I). If this should turn out not to be the case experimentally (see ||22||), 
the effect on the W mass may be correspondingly reduced. Here we just want to bring 
attention to the fact, that such a reduction cannot a priori be expected to depend propor- 
tionally on A, as it is presumably not the correlation strength, but the absolute number 
of correlated pairs of pions in the enhanced region of the Q distribution that is important 
for the size of the possible shift in the W mass. 

The comparison of hadronic decays of W + W~ — > qqqq and W + W~ — > qqlv channels 
gives a unique possibility to investigate the influence of Bose-Einstein correlations on var- 
ious properties of final state particles, such as multiplicity, transverse and longitudinal 
momentum spectra, resonance properties and reconstructed jet characteristics. It is pos- 
sible, that by taking proper care in the fitting procedures used, one may at the same time 
be able to use a large part of the hadronic W + W~ events for the W mass determination, 
and to study the interconnection effects in the relatively clean setting of e + e~ — > WW 
events, by restricting the study to the region of large, off-peak W masses where these 
effects are expected to be the largest. 

We would like to thank A. Olshevsky, S.A. Gogilidze, R. Lednicky, V.L. Lyuboshitz, 
T. Sjostrand and A. Tomaradze for helpful discussions. One of us (V.K.) would like to 
thank ALEPH collaboration for their kind hospitality at CERN, where a part of this work 
was done. 

Note added. Results of a similar study have been published recently p3| . The 



weight system used in |23| is a simplified version of our V4, where the computational 
difficulties were avoided by averaging over "clusters" of particles. The authors do not see 
any mass shift due to BE correlations at the level of their statistical precision. This is 
not inconsistent with the present calculation, as they have chosen to use R = 1 fm, which 
strongly reduces the number of interfering pairs. 
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